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The impact of filling up the Fermi sphere with the quarks, which dynamically generated their 
masses on the instanton liquid at finite temperature and baryonic/quark number density, is 
investigated. It is demonstrated, in particular, that the boundary of chiral symmetry restoration 
phase transition is shifted to the larger (about 100 MeV more) value of quark chemical potential 
compared to the magnitude inherent in the Nambu-Jona-Lasinio model. 
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Impressive results obtained in experimental study of ultrarelativistic ion collisions at RHIC 
(Brookhaven) and the experiments which are planned for the near future at ALICE LHC (CERN) 
PQ are standing in need of more accurate and precise theoretical predictions for possible signatures 
of new states of strongly interacting matter with fastly growing acuity. However, the theoretical 
advancement is much less appreciable especially in the latest years. For example, the predictions of 
various approaches for the behaviour of gluon condensate at finite T and non-zero values of bary- 
onic/quark chemical potential /i which is a key quantity for theoretical analysis are still inconsistent 
as before and at times simply conflicting. The possible changes appearing in the gluon sector at such 
conditions and usually described by varying the constants of multiquark interactions as the functions 
of T and \x in the Nambu-Jona-Lasinio (N JL) model [2] , need drawing almost inevitably the lattice 
numerical calculations of the gluon condensate [3] to be analysed. Practically to the same extent 
this remark is justified for the predictions of the chiral perturbation theory (CHPT) [3] and the 
\ QCD sum rules (SR) [5J. Both approaches have rather limited reliability for results of calculations 
^ ' around the critical parameter values. Actually, in Ref. [B] we have already tried to estimate the 
■ gluon condensate behaviour in hot and dense medium using the instanton liquid (IL) model [7j-[2] 
as an operative tool. In this case the screening impact of quarks filling the Fermi sphere up0 on the 
gluon condensate has been calculated for the massless quarks. 

In this paper we consider the influence of quarks with the finite masses on the gluon condensate. In 
the IL model the calculation of dynamical quark mass at zero temperature is grounded on making use 
the zero mode approximation |12j . However, even that calculation runs into rather serious technical 
difficulty (see also |13j ) while interpreting the loop quark diagrams at the chemical potential values 
exceeding the magnitude of dynamical quark mass \i > M q . We are treating this point here based on 
the NJL model and are not interested in the asymptotic large values of \i and T (see, for example, [Hj) 
and omit an analysis of the colour superconducting phase as well as the discussion of the difficulty 
in stabilizing the instanton ensemble which is rather often resulted in the speculations about the 
'realistic' structure of vaccum configurations. 



1 Following [TU] we take this effect as the dominating one at high temperature though there exist the other interesting 
possibilities |llj . 
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1 Approximating the vacuum conhgurations 
at finite T and \i 

Our purpose here is to find the practical and effective tool for evaluating the gluon condensate under 
extreme conditions. Obviously, such a task has been pending for rather long time and the prevailing 
number of scenarios to resolve it is grounded on the mean field approximation which supposes, 
actually, the simplified description of a system. It will be a guiding element of our approach while 
dealing with the instanton ensemble. For example, vacuum correlation function (A fl (x)A u (y)) of 
mean field description is transformed into the correlator which in the context of our approach leads 
to the mass generation for gluon field and, hence, to the colour screening factor. 

As well known in the IL model at zero values of p and T the superposition of (anti-)instantons 
in the singular gauge 

K{ x \l) = ~ ^ ab Vb^u a v (y) , a u (y) = — — f , y = x-z, p, v = 1, 2, 3, 4 , (1) 

M g y +p y 

(where p is a pseudo-particle size, u is a matrix of its colour orientation and z is its center coordinate) 
is considered as the ground vacuum field saturating the QCD generating functional (dealing with 
anti-instantons one has to change the't Hooft symbol fj — > rj). The QCD generating functional is 
evaluated to be as 

N oo 1 N 



OO 1 JV „ OO 1 iV . 

Y = E II d o(Pi) UM = E 4 II ** e- E( " ] , (2) 

7V=l iV - i=l J N=l iV ' i=l J 

E(j) = (3 U int {i) -J2 lndo(pi) , i = (z,p,u), 



here 



d ( P ) = j (3 2N < e-M (3) 

is the distribution function over the size of individual instanton (dilute instanton gas approximation) 
[7], d'ji = d A Zi duJi dpi is the integration element, 

M = K = -* Hc^Ap) 
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is the single instanton action (A = A^j^ = 0.92Apy) with the constant depending on the 
renormaWion SC heme C N « ^£t^% with alTOther parameter \ = HHszMl. 

The auxiliary coefficients (3 = —b ln(Ap) (3 in the exponent of Eq. (J2J) are fixed at the characteristic 
scale p (pseudo-particle average size). Assuming the topologically neutral instanton liquid we do not 
differ the instantons and anti-instantons and denotes (when used) the total number of pseudo- 
particles which occupy the volume V. 

Taking into account the interaction of instantons with vacuum fluctuations is effectively presented 
by appearance of the screening factor in the distribution (jBJ) 

d{p) = 1 ^ 6 -/3(p)-Cp 2 ( 4) 
P 

where the magnitude of screening coefficient ( is dependent on the choice of superposition ansatz. 
For the pseudo-particles in the singular gauge the interaction term taken in the pair interaction 
approximation is [9] 

rf^i du 2 dz l dz 2 U int (71,72) = V £ 2 p\ p\ , 
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with the constant £ 2 = 



Nr. 



-. It is interesting to notice here that the configurations used in 



c 

the valley method [15] result in the significantly smaller value (about one order) of coefficient £ |16j . 
Besides, the screening factor can be steadily extracted from the lattice data as ~ 0.22 fm [TTj with 
the configuration cooling procedure. The corresponding configurations are reasonably well fitted by 
the instanton ensemble as was shown [IB] , although the analysis of optimal instanton configurations 
in the mean field approximation is worthy of special study and will be done in the separate paper 

The convexity property of exponential function allows us to estimate the partial contribution into 
the generating functional Eq. (T5]) at each value of N by the following approximating form 



Y > Y 

1 / 1 r. 



approx 



Y 1 exp(-(E~E 1 )) , 



which can be presented [9] as 

Y a pp rox 6 , X N 



V 



+ 1 [ln(n/A 4 ) - 1] - JVln 



(5) 



(6) 



where n = N/V, v — (b — 4)/2. Then the respective parameters of IL are defined by the maximum 
in n of generating functional with the interrelation of instanton average size and its density taken 
into account 

K 2 np 



v 

7i2 



Now calculating the maximum of X in n we have to resolve the following equation 



+ 1 ln(n/A 4 )+ln 



2N c d(3 v dp 
+ n — n —^.—r- 

(3 an 2p an 



(7) 



(8) 



Owing to the relation (J7D we have 



ld(5 Idn 4 
p dp n dp p 



From the other side 



d/3 



b dp _ dp 



dp 



P' ~dp = ~dp' -^- ewr ^ m S ^ ne derivative of P in the density as ^ 



ip 



^/^p, we come to the system of equations 



dp 
dn 



1 bp 
n 4p-b 



dp_dp 
dn dn 



(9) 



Finally resolving the system of transcendental equations we can determine the equilibrium IL pa- 
rameters. 

At the finite temperature the configuration saturating the generating functional is changed by 
the superposition of (anti-)colorons [20], [21] which are the periodical in the Euclidean 'time' (with 
the period of T _1 ) solutions of the Yang-Mills equations [22J i.e. 



<$>(x,T) 



1 + 



np 2 T 



sinh(27rrT) 



(10) 



r cosh(27rrT) — cos(27ttT) 



Here r = \x — z\ defines the distance from the coloron center z in three-dimensional space, r = x^ — 
is the 'time' interval. It can be easily seen that the solution is transformed into the (anti-)instanton 
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one in the singular gauge at temperature going to zero. Clearly, the distribution function over the 
coloron size [10], [23] is also changed 



d(p; fi, T) = d{p) p 2 , rffr, T) = 2 vr 2 



Nr. 



A' 



6 /=i 



The first term of the screening factor describes the one-loop gluon contribution into the effective 
action and the second term generated by quark contribution in one-loop approximation can be 
exactly calculated and is free of the 'bad' singularities [24] . The 'time' component of polarization 
tensor generated by quark of fixed colour has the form 



n{ 4 (fc 4 ,w) 



k 2 



h 2 uj 2 



•pijj 



dp p 2 



n. 



1 + 



l::-k 2 _ {k 2 + 2pu) 2 + 4e 2 p k 2 



8pk 



In 



(k 2 -2 P Lu) 2 + 4e 2 p k 2 



arctan 



8pu> Epk^ 



Aeiki - ApW + k 



here u 
(e^ + 1) 



\k\, k 2 = u) 2 + k%, e p = (M 2 + p 2 ) 1 / 2 where M q is the quark mass 



, Tip 



n„ + nj. 



n. 



n 



e t 



1) _1 (n~,n p are the densities of anti-quarks and quarks, respectively). 
When summed up over all the components the polarization tensor can be presented in the following 
form 



n^(fe 



4- 



7T 



dp p 



n, 



1 



2M 2 



8pk 



k\ (k 2 + 2pu) 2 + 4e 2 k 2 
In 



(k 2 -2pcu) 2 + Ae 2 p k 2 



'12) 



It is clear when the zero-component k± = the dominant contribution into the gluon mass at small 
values of u comes from the first term (a unit) and the space components are negligible. In particular, 
at ui = it will be 



n^,T) = n{ 4 (o,o) 



7T 



dp p 



(13) 



{p 2 -M 2 fl 2 p M 2 ^p + (p 2 -M 2 y/ 2 



In 



M„ 



. In order to calculate 



Then at T = we have II / (u, 0) 

the IL equilibrium parameters as the functions of p and T one has to minimize the approximating 
functional © making the substitutions of ((7j) and © for 



= = V 2 + K 2 np 2 
p 2 



(14) 



and 



ndf3 b 
]3dn~ 4(3 -b + 



(15) 



correspondingly. 



2 Quark mass generation in stochastic field 

It is anticipated that in the IL model the quarks are considered to 'live' (and to be influenced) in 
the stochastic (anti-)instanton ensemble which is defined by the following generating functional 

Z = jD^]D^](eI dxC ") A , C q = 4>{x) [id x + j^gA{x ]lk ^{x) , (16) 
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where the averaging over (anti-)instanton ensemble is implied. The consistency requirement for 
effective Lagrangian in the Hartree approximation results in the equation for the quark Green function 
which reads as 



M(p) = -y i £ / ft tJ^I (2vr) 4 5 4 (t ft) Tr (g A( Ql ) S(p -qj.-.g A(q n 



(17) 



Being summed up the right hand side of Eq. (ITTj) can be presented in the compact form as [26] 



N 



M(p) 



N„y 



i=i 



TrJ2 (p| 



-i 



-i 



[P> 



18) 



(in such a form the averaging over the pseudo-particle location z and calculation of colour trace is 
meant). Analyzing the solution in the form 



S(p) 



p — iM{p) 



(19) 



where M(p) denotes the quark mass, one can calculate the highest term of expansion in the IL density 
(presented by the zero quark mode $(p) in the instanton field) as 



M(p) ~ n 



1/2 



p 2 $ 2 (p) 



1/2 



At finite quark chemical potential the derivative id in Eq. (fT6l) should be substituted for id — ifi 
where fx = ^74. Then quark Green function ( fl9l) develops the following form 



where 



M(p; li) ~ n 



S(p]jj) 



1/2 



1 



p + ifl — iM(p; n) 



(p + tLi) 2 $ 2 (p;/i) 



(20) 



dp 



(p + iLi) 2 $ 2 (p;/i) $ 2 (p;/x) 



1/2 



and the overt expression of the zero mode could be found, for example, in |27J. With the chemical 
potential increasing and reaching the values of dynamical quark mass order (ll ~ M q ) the quark 
mass magnitude M(p; ll) begins to increase as a power. The similar situation with the dynamical 
mass increase takes place for the approach in which the unperturbated quark Green function So is 
approximated by the zero modes [28p . Such a behaviour is non-physical and contradicts to the 
intuitive expectations. It seems, the situation could be improved by taking into account the non- 
zero mode contributions but very complicated analytical structure of the corresponding expressions 
makes this calculation practicaly hopeless. Thus, the question about the estimate of non-zero mode 
contribution is still vague (see, however, [29]). The proposition to treat poles as in Ref. [13] leads, 
unfortunately, to unphysically small values of the chemical potential of chiral symmetry restoration 
phase transition. It is interesting to notice here that the zero mode approximation is quite reliable 



2 Determining the saddle point parameter in this case one encounters the problem of calculating a loop integral in 
which the pole of Eq. (pUf appears on the real axis. The treatment of that integral as a principle value gives its real 
part only which is not enough, of course. However, this problem is softened by the situation that the pole appearance 
on the real axis occurs in the local vicinity of transition point into the colour superconducting phase and, therefore, 
more precise definition of the saddle point parameter looks superfluous. However, the general problem of calculating 
the loop integrals is entirely hot and actual in the context of analyzing the chiral symmetry restoration. (Prof. T. 
Hatsuda has drawn our attention to this aspect of the problem and the authors are very grateful to him for that.) 
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even at the finite temperature if one confines oneself to work with the chemical potential values not 
larger than the dynamical quark mass fi ~ 300 MeV. The quark condensate estimates are quite 
suitable in this case even if one makes use simply non-coloron zero mode. These results put forward 
the task of searching the effective approximations for the equations of type (ITTj) . The significant 
progress in studying the systems of quarks at finite temperature and chemical potential has been 
reached in the framework of NJL approach. (Let us remember here that at finite temperature the 
integration over the fourth component in Eq. fll6p should be performed in the interval from zero 
till T~ l and gluon fields obey the periodic boundary conditions whereas the fermion fields obey the 
anti-periodic ones.) 

As it is difficult to handle (anti-)instanton ensemble directly we are going to retain some essential 
features of (anti-)instanton configuration contribution and approximate it with the simplest form. 
Actually, we suppose the existence of superposition of stochastic randomly oriented color gluon fields 
in the Euclidean space. These fields have the (^-function form with their randomly distributed centers 
z, i.e. 

A^x) = U ] r a U (2vr) 4 5 {4 \x - z) , A„{p) = U ] r a U (2vr) 4 e ipz = e ipz . (21) 

It is clear if one considers one-particle correlations only (just what is done for the pseudo-particle 
ensemble) the simplest non-trivial correlation function (A(xi)A(x2)) z ,u will lead to the point-like 
interaction of quarks 

(A(x 1 )A(x 2 )) z , u ^5^(x 1 -x 2 ) 

which is specific for NJL. In further analysis we do not need to know the concrete form of stochastic 
factor A and do not specify it here. Searching the solution of Eq. ( JT71) in the form M(p) = M we 
introduce the auxiliary function ip{q.,p) the following way 



M = 4^ f T^Ti Tr (A(q - p) 1 — j>(q,p)) M ,u 

N c J (2tt) q + in, — iM 

(22) 

f dq 1 

tp(p, p') = A(p -p')+ / ——4 A(q - p) — — — ip(q, p') , 

J {2tt) q + ijj, — iM 

(of course, we imply non-zero quark chemical potential). Presenting the solution for ip(q,p) in the 
form tjj(q,p) = ip e^ g ~ pS>z we are able to obtain the following equation to determine the function ip of 
our interest [f| 

ij = A + A [ -^j - — ^ — - v • (23) 
y J (2tt) 4 q + ijl-iM y v ' 

As in the NJL model Eq. (1231) requires the regularization. Here we are using the conventional 
procedure of three-dimensional momentum regularization [2] which allows us to obtain 

/= / 7^14 rr - l . M =iC lA + iDM , (24) 
J (2tt) q + ifji — iM 

6{y - M) {fi 2 - M 2 f' 2 

(2tt) 2 3 



3 Another utmost regime where the correlation length is supposed to be infinitely large (A(x)A(y)) — A 2 is also 
very interesting. It was analyzed in the Keldysh model [30] and the exact solution was found out. The complete 
summation of series for the quark Green function results in the expression as 

1 f°° -t 2 1 



(2tt) 1 /2 p + fi-tA' 

(it is given here in the Minkowski space). Apparently, it has no poles similar to the non-relativistic Green function as 
well as the 'analytical' model of confinement 3 lj . 
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D = < 



■57T' 



8vr 2 



A+(A 2 + M 2 ) 1 / 2 
M 



A (A 2 + M 2 ) 1 / 2 -M 2 ln 
A (A 2 + M 2 ) 1 / 2 - fi (fi 2 - M 2 )V2 _ m 2 In 



, fi<M 
A + (A 2 + M 2 ) 1 / 2 



n > M 



A denotes here the cut-off value of three dimensional momentum in the I integral. Finally, we have 
for the solution of Eq. (1251) the following result 



B + C 2 A 2 
B(1 + D 2 A 2 M 2 ) 



(A + iDA 2 M) 



iC 



B(1 + D 2 A 2 M 2 ' 



(A + iDA 2 M) 74 (A + iDA 2 M) , (25) 



where B = 1 — 2iCA i — C 2 A 2 + D 2 A 2 M 2 . Using this solution in Eq. (1221) . averaging over the colour 
orientation and holding the highest terms of the N c expansion we come to the mass gap equation 



M = An 



DA 2 M 



1 + D 2 A 2 M 2 - C 2 A 2 ' 



(26) 



It enables to formulate the condition which signals the breakdown of chiral symmetry (the generation 
of quark mass) if such a constraint is obeyed 



A 2 > 



1 



An D — D 2 M 2 + C 2 ' 



(27) 



If one neglects the contributions proportional to A 2 in the denominator of Eq. ( 1261) the gap equation 
of the NJL model with the coupling constant G of four-fermion interaction is exactly reproduced 



G = An A 2 



A 2 



57T 



(28) 



In order to receive the qualitative estimates it is worthwhile using the characteristic cutoff parameter 
of the NJL phenomenology A ~ 600 MeV. In our estimates we rely on the constraints D 2 M 2 A 2 1, 
C 2 A 2 < 1 only. 
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Figure 1: The effective potential W as the function of quark mass for different values of chemical 
potential fj, = (lower curve), fi = 100, /i = 200, [i = 300, fi = 400 fi = 450, fi = 500 MeV. 
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The estimate of A 2 for instanton ensemble could be obtained from the corresponding correlation 
function [32] 

(A{x - z) A{y - z)) Z)U = F{x - y) . 

Actually, we have F(0) = jjlr- — t p 2 , which means A 2 . ~ \ and then we receive for the Fourier 

iV c i p 

component A 2 ~ p 6 . Taking into account that p ~ A -1 and D ~ A 2 (see Eq. f l24|) ) we obtain 
D 2 A 2 M 2 ~ '^ ie s ^ anc ^ arc ^ parameter values of the NJL model provide us with the small 

magnitude of this factor and it could be neglected. Surely, it is a fairly serious argument in favour 
of using the developed approach. 



3 Approximation of NJL and IL model 



Our above analysis demonstrates that approximating the instanton correlator with the delta- function 
form and using the regularized NJL model at the same time is the fully compatible procedure. 
Moreover, Eq. f[2"H]) provides us with the possibility to consider the interrelation of gluon and quark 
sectors. It relates the constant G and the IL parameters such as the IL density n and average 
potential A and G is related to the dynamical quark mass M = M(0) which defines the screening 
effect. 

Let us now remind that the generating functional of the NJL model has the following form [2] 



Q = Gn(i z 



N f N r 



7l~ 



p 2 e p dp 



N f N r 



-T 



7T 









In ^ 


lo 





p 



1 + e T 



+ In 



e p — p 
1 + e T 



(29) 
dp, 



where e p 
we have 



(M 2 +p 2 ) 1//2 , M = m — 2Gq a, m is the current quark mass and for the quark condensate 



-M 



A P 2 
6 V 



n 



dp , 



a 



(30) 



The quark mass is defined by calculating the minimum of Q as the function of M (or quark condensate 
a) and the coupling constant Go together with the cutoff parameter A is fixed phenomenologically (by 
fitting the experimental data). We suppose to take this quantity as an estimate of quark determinant 
while quarks are in the stochastic field of (anti-)instantons and modify the determinant aiming to 
include the interrelation of quark and gluon sectors. In this way we use (instead of G) in Eq. (129]) 
G — » Go where n , Go are the IL density and the constant of fourquark interaction at zero 
temperature and zero chemical potential. In full analogy with the IL model it is easy to understand 
the parameter A should generate the factor of the p 3 type and the substitution of the coupling 

constant for 



n P 6 Go looks quite natural. On the other hand as the simplest option we could take 
n oPo 



the cutoff parameter in the quark sector A unchanging as the parameter A 
quarks in stochastic instanton ensemble we should change Q in Eq. ( 1291 for 



A" 



Thus, for the 



n 



Q = —Go^ 

no 



NfN c 



TV 



p 2 e p dp- 



N f N c 

7T 2 



T 



V 



e 



In 1 + e 



P + P^ 



+ In 1 + e 



~'P P 



dp. (31) 



Apparently, the vacuum parameters of the IL model and the NJL one should not change. In order 
to realize that one should make the corresponding subtractions just to retain the effect caused by 
the quarks filling the Fermi sphere up because the interrelation of vacuum (at zero T and p) quark 
and gluon fields has been already discounted effectively in the running coupling constant and by 
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Q Q I I I I I 

' 200 400 600 800 1000 

M (MeV) 

Figure 2: The IL density as the function of chemical potential for different values of temperature 
(the IL density is decreasing with the temperature increasing) T = 0, T = 50, T = 100, T = 150, 
T = 200, T = 250, T = 300 MeV. 

tuning the NJL model parameters. The (anti-)instanton ensemble and quark field are described by 
the product of functionls Y approx and Z. Therefore, for the effective potential we have 

w = x+n . 



0.7 r~i — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — r 




* 100 200 300 

T (MeV) 

Figure 3: The IL density as the function of the temperature T. The lower curve corresponds to the 
calculation with massless quarks. 

The equilibriumn IL parameters are defined by the effective potential minimum on the IL density, 
i.e. by ^p- = 0. However, as was declared the following subtraction 



should be done. Similar operation should be executed at determining the quark mass = 0, i.e 



dX dr] 2 dX drf 



drf dM drf dM 



\lx=0,T=0 



+ 



on 

dM 







(33) 



Two first terms of Eq. (133]) are the result of the fact that the overt dependence on the quark mass 
in the contribution into the effective potential is available in the screening factor rf only. Strictly 
speaking one should integrate till the momentum order of A in Eq. ffT3l , too. However, such an 
amplification is superfluous as the detailed analysis shows. In practice, obviously, it is simpler not 
to resolve the transcendental equation (which has two branches, at least) but calculate simply the 
minimum of effective potential in M. It is easy to understand that for the concrete form of our quark 
effective potential Q the equation for determining the equilibrium IL parameters coincides with the 
vacuum one because the direct dependence on the IL density n is present in the first term of Eq. (I3"T|) 
only. Just because of that reason two last terms of Eq. (!32l) are canceled. The subtraction in the 



second equation of (133]) should not be performed as 



dX 



n 



o 

fj- and the function rf = at i± < M. 



drf 

Thus, the equilibrium IL parameters are defined by the same scheme as before [6] and minimum of 
the generating functional W in M fixes the dynamical quark mass. 

Here we are using the following set of the NJL model parameters [2] (T. Hatsuda, T. Kunihiro). 
We take for the current mass of u and d quarks the same value m = 5.5 MeV, for the cutoff parameter 
as A = 631 MeV and for the ratio of the coupling constant to its critical value as a = Gq/G c = 

1.33/Nf, ( G c = ] • Such a set of parameters results in the following values of the 7r-meson mass 



N C A\ 

139 MeV and the constant of pion decay F n = 93 MeV. 
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Figure 4: The dynamical quark mass as the function of chemical potential for different temperature 
values T = (upper line), T = 50, T = 100, T = 150, T = 200 MeV. 



The Fig. [T] shows the effective potential W as the quark mass function for the different values of 
chemical potential fi. The lower curve corresponds to zero value of chemical potential and is in full 
coincidence with the respective curve of the NJL model. With the chemical potential increasing the 
process of filling up the Fermi sphere starts and it is easy to see that the effect of pseudo-particle 
field screening by the quarks of small masses occurs dominating (see Eq. ( TlBl) and the next one). 
The screening effect leads to diminishing the absolute value of gluon condensate and, hence, the 
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effective potential of (anti-)instanton and quark system is increasing. It is distinctly visible in Fig. 
[T] in the region of small quark masses. With the quark mass increasing, the impact of the filling 
up process is amplified. Starting on some value of small quark mass the saturation regime manifest 
itself and the gluon condensate is suppressed. The plateau is formed and it is well seen in Fig. CD At 
the larger values of chemical potential the chiral symmetry restoration starts and the corresponding 
curve is shown with the dashed curve in Fig. [TJ More detailed analysis allows us to conclude the 
simplifications made do not depreciate the qualitative picture of screening. 
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fi (MeV) 

Figure 5: The quark matter density as the function of chemical potential. The last right curve 
corresponds to the zero temperature. The next left curve corresponds to the temperature 50 MeV, 
then 100 MeV, and eventually T=150 MeV. rio q = 0.062 fm -3 is the normal density of quark matter 
calculated from the normal baryon matter density ns = 0.45 fm -3 . 

The IL density as the function of chemical potential at the various values of temperature (with 
the temperature increasing the IL density is decreasing) is plotted in the Fig. [21 The dashed curve 
which was obtained by us in [6] (the similar mechanism of screening was discussed also in [TO], |20j ) 
corresponds to the calculation with the massless quarks. The mechanism of forming the observable 
plateau is quite understandable. Until the chiral symmetry restoration does not take place the quarks 
are ineffecive in the gluon field screening and the gluon condensate practically does not change. 

Fig. El presents the dependence of IL density on the temperature. As it was expected the density 
is larger for the massive quarks than for the massless ones in the region of low temperature (below 200 
MeV). This result agrees qualitatively with the observation done in [33]. Some lattice calculations 
support this scenario of screening. For example, in [31] it was proven that the Debye screening mass 
behaves as m e i ~ gT and depends on the quark flavours (see Eq. (jTTT) ). Exponential suppression of 
gluon field with increasing temperature was also found out in the lattice measurements of correlation 
functions dealing with the cooled configurations [35]. In fact, gazing into the detailed analysis of the 
problem under consideration we collected a lot of reasons to have the topological solution with the 
suppressed chromoelectrical component instead of the (anti-)coloron one to construct more realistic 
approach. In particular, it was noticed in [SB] the coloron solution mainly does not fit the lattice 
data. However, we understand such global pretension (as a disproving conclusion) is rather naive 
because our result here shows the coloron solution is quite practical for estimating the screening 
effect. 

Analyzing the quark sector we calculated the behaviour of dynamical quark mass as the function 
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Figure 6: The ratio of quark condensate value to its magnitude at zero temperature and zero chemical 
potential as the function of temperature at the various values of chemical potential /z = (last right 
curve corresponds to the zero value of chemical potential), \i = 110, \i = 200, /x = 300, /i = 400 MeV 
(last left curve). 

1 200 i i 
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Figure 7: The masses of 7T- and a-mesons as the functions of chemical potential at different values 
of temperature T = 0, T = 50, T = 100, T = 150 MeV. The dashed lines correspond to the a-meson 
and the upper dashed line corresponds to the zero temperature. Last right solid line shows the 
7r- meson mass behaviour at zero temperature. 

of chemical potential and plotted it in Fig. H] for various values of temperature. The upper curve 
corresponds to zero temperature and the quark mass behaviour along this curve coincides with the 
NJL model up to the chemical potential value of /i ~ 300 MeV. With a further increase of /i the quark 
mass increases. It is quite understandable qualitatively if one looks at Fig. [TJ At the commencement 
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of chemical potential increase the screening effect does not produce any noticeable impact on the 
minimum of effective potential W. In spite of effective potential increase in the region of small quark 
mass the threshold value of chemical potential should be reached at which the forming plateau begins 
to expel the effective potential minimum to the larger mass values. The size of region in which the 
quark mass increase takes place is comparable with the quark mass (order of 100 MeV) and is of 
interest, in particular, for investigating the equation of the state of strongly interacting matter. In 
the version of the NJL model with the parameter choice suggested by T. Hatsuda and T. Kunihiro 
[2] the chiral symmetry restoration occurs at quite low density (the order of normal nuclear matter 
density). In our approach the significant decrease of quark mass is shifted (drags on) to the region of 
large chemical potential values approximately 100 MeV larger, which agrees entirely with an intuitive 
expectation. Fig. [5] is devoted to analyzing the quark matter density as the function of chemical 
potential. It seems the shift of chiral symmetry restoration phase transition to the region of larger 
chemical potential values (order of 400 MeV), could generate an essential increase of quark matter 
density. However, Fig. [5] demonstrates the change in this interval is inessential (the increase of quark 
mass n q ~ (/i 2 — M 2 ) 3 / 2 provides the compensation), and in actual fact we have to deal with the 
same vacuum quarks as at // < 300 MeV. 

The quark condensate normalized to its value at zero temperature and zero chemical potential 
is depicted in Fig. E] as the function of temperature for different values of chemical potential. The 
behaviours shown are in full agreement with the predictions of the other models. Finally, two last 
Figs give more information on the masses of tt- and a-mesons which are also calculated in the NJL 
model (see, for example, M.K. Volkov, A.E. Radzhabov [2]). The 7r-meson mass is given by 

*J = rf-55j?. (34) 

where g^ qq = is the renormalized constant of meson field interaction including the following 
auxiliary integral 

/2 = 8^ Jo 4 {1 - n *> ~ np)dp - 
The mass of cr-meson is defined by the mass of 7r-meson and dynamical quark mass as 

M 2 = M 2 + 4M 2 . (35) 
Pion decay constant which is a key element of model tuning is defined as F~ = Fig. [7J 

iJnqq 

presents the masses of it- and cr-mesons as the functions of chemical potential for various values of 
temperature. The dashed curves correspond to a-meson. The upper dashed line shows behaviour at 
zero temperature and the solid lower line corresponds to zero temperature behaviour of the 7r-meson 
mass. The interval in which the a- and 7r-meson masses become identical defines the parameters (on 
the fi-T plot) corresponding to the chiral symmetry restoration. It is clear from Fig. [7J that such a 
restoration at zero temperature occurs around [i ~ 460 MeV, at T = 50 MeV around \x ~ 410 NeV, 
at T = 100 MeV around \i ~ 350 MeV and at T = 150 MeV around // ~ 220 MeV. Besides, this plot 
allows us to fix the line m CT = 2m n on which the strong decay channel of a-meson is close. At low 
temperature (T < 100 MeV) the 7r-meson mass undergoes a significant change in the region of chiral 
symmetry restoration only and the estimate that the line m a = 2m n approximates chiral symmetry 
restoration curve (m a = m n ) looks rather practical. At T = 100 MeV the chemical potential for the 
line m a = 2m n is about fi ~ 320 MeV and at T = 150 we have fi ~ MeV. The details of this line 
behaviour could be quite indicative for searching the mixed phase in relativistic heavy ion collisions 

mi- 
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4 Conclusion 



In the present paper we investigated the effect of gluon condensate screening with the massive 
quarks filling up the Fermi sphere. We developed the approach based on the NJL model highlights 
which allows us to get informative qualitative estimates. In particular, we argue that one of the 
manifestations of filling up the Fermi sphere could be an increase of the quark mass and, hence, 
the shift of chiral symmetry restoration phase transition to the larger values (about ~ 100 MeV) of 
quark chemical potential. Another instructive result obtained implies that the gluon condensate does 
not die out completely in the parameter region characteristic for this phase transition even at the 
most advantegeous regime of vacuum gluon field screening. The lattice measurements of the same 
quantity confirm such a conclusion allowing us to predict that gluon condensate are surviving even 
in the region of parameters essentially beyond the values admissible for our approximation. 
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Efimov, Yu.B. Ivanov, Yu.L. Kalinovsky, A.E. Kuraev, A.E. Radzhabov, V.V. Skokov, V.D. Toneev, 
V.L. Yudichev for numerous fruitful discussions. Financial support within the Grant INTAS 2004-398 
is greatly acknowledged. 
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